Spin dynamics of the model 2D quantum antiferromagnet CFTD 
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The magnetic excitation spectrum in the two-dimensional (2D) S — 1/2 Heisenberg antiferro- 
magnet copper deuteroformate tetradeuterate (CFTD) has been measured for temperatures up to 
T ~ J/2, where J = 6.31 ± 0.02 meV is the 2D exchange coupling. For T< J, a dispersion of the 
zone boundary energy is observed, which is attributed to a wavevector dependent quantum renor- 
malization. At higher temperatures, spin-wave-like excitations persist, but are found to broaden 
and soften. By combining our data with numerical calculations, and with existing theoretical work, 
a consistent description of the behaviour of the model system is found over the whole temperature 
interval investigated. 
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While at the atomic level, magnetism is a quantum 
phenomenon, the collective behavior of magnets can 
be largely understood using classical concepts. This 
can even be true for antifcrromagnets with low spin 
and spatial dimensionality, where quantum fluctuations 
are sizable. A famous example is the 2D quantum 
(S = 1/2) Heisenberg antiferromagnet on a square lattice 
(2DQHAFSL) with nearest-neighbour interactions. Con- 
siderable effort has been devoted to this particular model 
because it describes the parent compounds of the high-T c 
superconducting cuprates, and also because it once was 
thought to have a spin fluid rather than a Neel ground 
state. 

By now, it is well established that at zero temperature 
the 2DQHAFSL displays long-range order, albeit with a 
staggered moment reduced by quantum fluctuations [0. 
In addition, harmonic spin-wave (SW) theory, based on 
a classical image of the spins as coupled precessing tops, 
gives an excellent account of the spin dynamics up to in- 
termediate frequencies at T=0. At finite temperatures, 
long-range order is destroyed by thermal fluctuations, 
and the system possesses short-range order only, char- 
acterized by a temperature-dependent correlation length 
£(T). A ombination of low temperature renormalised 
classical theory 0], intermediate temperature quantum 
Monte Carlo (QMC) MM and high-temperature meth- 
ods |M accounts for the experimental data for £(T), 
covering the range J/5 < T < J pj-|To|]. Thus, a co- 
herent picture exists for the thermodynamic properties 
at all temperatures for S=l/2, as well as for higher val- 
ues of the spin [ pd| . The remaining questions about the 
2DQHAFSL concern the intermediate and low frequency 



spin dynamics at nonzero T, and the high frequency spin 
dynamics at all T. In both cases one could well expect 
to see more severe quantum effects because they should 
be more sensitive than the static properties to the non- 
linearities of Heisenberg's equations for a low spin sys- 
tem. The present paper describes the first experiment 
to confront all of these questions directly over the full 
energy scale set by J. 

Though experiments on the dynamics of the 
2DQHAFSL are relatively scarce, considerable theoret- 
ical work exists. Time dependent information is con- 
tained in the dynamical structure factor S(Q,uj) = 
f duje- iwt J2r eiQr ( S o(°) S r(t)), which is a function of 
wavevector Q and energy hw, and is measured directly 
in neutron scattering experiments. The T — proper- 
ties can be described by classical SW theory, but with 
quantities such as the SW velocity v s = Z c \J~2Ja, spin 
stiffness p s — ZpJ/4, and susceptibility XX = Z X /{^,J) 
renormalized by quantum corrections. Consistent values 
for the renormalization constants Z c = 1.18, Z p = 0.72 
and Z x = 0.51 have been obtained using the SW expan- 
sion to order 1/(2S) 2 Q, and series expansion from the 
Ising-limit jfjj. However, these two approaches to the 
T = spin dynamics disagree in detail, in that the SW 
expansion predicts a constant energy at the zone bound- 
ary (ZB), whereas the Ising-limit expansion predicts a 
~7% dispersion. 

On warming, scattering from thermally excited 
magnons, and a variation of the order parameter on a 
length scale set by £, will limit the lifetime of the exci- 
tations. Correspondingly, the inverse magnon lifetime, 
r, will increase, while the dispersion softens to lower 
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energies. Various theoretical studies ]13|-|20| and QMC 
calculations in the limited temperature range 0.35J < 
T < 0.5 J [pl| have addressed the dynamic structure fac- 
tor of the 2DQHAFSL at finite T. Most approaches 
agree on the existence of well-defined excitations, but 
few experimental data are available to test the vari- 
ous quantitative predictions. Much of the experimen- 
tal work on the 2DQHAFSL has focussed on La2Cu04 
§ and Sr 2 CuCl 2 2 [§j2|, both of which contain the 
Cu0 2 planes that are the building blocks of the cuprates. 
These materials have the drawback that their large cou- 
pling constant J ~ 1500 K makes studies on temper- 
ature and energy scales comparable to J technically 
challenging. Another realization of the 2DQHAFSL is 
Cu(DCOO) 2 4D 2 (CFTD), where the 5 = 1/2 Cu mo- 
ments are coupled through formate groups, leading to an 
exchange energy (J = 73.2 K) which is more amenable 
to experiments. The correlation length £(T) in CFTD 
has recently been measured up to temperatures of T ~ J 
H|o|, and found to be in good agreement with theory 
and computations. 

Below 236 K, CFTD has the P2\jn space group with 
lattice parameters a = 8.113 A, b = 8.119 A, c = 12.45 A, 
and a monoclinic angle (3 — 101.28°. The ab plane con- 
tains face centred 5 = 1/2 Cu 2+ ions, forming an al- 
most square lattice. The SW dispersion along a* and 
b* measured by neutron scattering is well described by 
an isotropic nearest-neighbour coupling J = 6.31 meV, 
and a small anisotropy induced gap of 0.38 meV at the 
zone centre |23 . The inter-layer coupling is estimated to 
be less than 10~ 4 J, while a Dzyaloshinskii-Moriya term 
Jd = 0.46 meV has been inferred from magnetization 
measurements p3. Below Tjv = 16.54 K the system or- 
ders three-dimensionally due to the inter-layer coupling 
& 

Our neutron inelastic scattering experiments using the 
direct geometry time-of-fiight (TOF) spectrometer HET 
at ISIS, UK. The incident energy Ei = 25 meV gave 
an energy resolution of FWHM 1.64 meV at zero energy 
transfer. A 3.71 g single crystal of CFTD was aligned 
with the 2D planes normal to the incident beam. In the 
geometry used the horizontal and vertical detector banks 
measured scattering along the [l,l]-type direction, while 
the diagonal banks collected data along the [1,0] direc- 
tions. The measured intensities were corrected for de- 
tector efficiency, and converted to absolute units by nor- 
malising to the incoherent scattering from a vanadium 
standard. 

Data were collected at several temperatures between 
8 K and 150 K with typically 24 hours of counting per 
temperature. The data at 8 K (within the 3D ordered 
phase) provide a precise determination of the spin Hamil- 
tonian, and allow characterization of the phonon back- 
ground. For example, the band at 20 meV is probably 
due to the localised motion of the crystal bound wa- 
ter, which freezes in an anti- ferroelectric transition at 



236 K ~ 20.3 meV |2q] , while the extra scattering around 
7 meV is due to acoustic phonons emerging from the 
(1,0,1) type reciprocal lattice points, see Fig. 0(a). 




FIG. 1. Inelastic TOF data in CFTD with Q 2D in units of 
(it, 7r). (a) Raw data at 8 K. (b) Representative cuts along en- 
ergy (bl) and Q (b2) as indicated by the windows in panel (a). 
The vertical scales cover to 100 mbarn/steradian/meV/spin, 
and the solid lines are fits to the resolution convoluted 
cross-section described in the text. (c)-(d) Background 
subtracted data at 16 K and 36.2 K respectively. The 
pseudo-colour scale (from blue to dark red) range from to 
50 mbarn/steradian/meV/spin. 

At 8 K, the magnetic signal is concentrated along a 
sharp SW dispersion curve. The data were analysed by 
taking cuts along Q and u>, as illustrated in Fig. |l|(b), 
allowing the phonon contribution to be dealt with as a 
local background. Each cut was fitted with a parameter- 
ized model for the scattering obtained by convolving the 
full instrumental resolution with the linear SW theory 

form S(Q, uj) = 



1-7Q 



1+7Q 



S(u - ujq); ujq = 2JJ1 



7q = ^(cosQ x + cos Q y ), multiplied by the magnetic 
form factor of free Cu 2+ and the Bose population factor. 
Since the parameterization is only applied locally (the 
parameters A and J were allowed to vary for each cut), 
it imposes no constraints on the global form of S(Q,lu) 
or on the extracted dispersion, shown in Fig. 2. It is 
evident that along [1,1] the dispersion is well described 
by a uniform renormalization of the linear SW theory re- 
sult. A fit gave J = J/Z c = 6.31 ± 0.02 meV, in good 
agreement with previous neutron scattering [E3| studies 
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and the value derived from high-temperature suscepti- 
bility |24|]. The present data can be combined with the 
susceptibility data to give a value of the renormalization 
factor Z c = 1.21 ± 0.05, within error the same as the 
theoretical value 1.18 lElO. 




(71/2,71/2) (7C,7C) (271,0) (71,0) 

FIG. 2. The SW dispersion in CFTD at 8 K. The solid line 
is a fit of the Q||[l, 1] part to nearest neighbour SW theory, 
giving J = 6.31 ± 0.02 meV. The dashed line is the result of 
an expansion from the Ising-limit J5) using J = 6.31 meV. 

Along [1,0] there is a pronounced dip in energy around 
(71-, 0), indicating a dispersion of 6 ± 1% along the zone 
boundary (ZB) from (7r/2,7r/2) to (vr,0). Such a ZB 
dispersion has also been observed in Sr 2 Cu304Cl2 [ffiH , 
where it was explained in terms of the Ising-limit ex- 
pansion |J. Indeed, the Ising-limit expansion accounts 
perfectly for our data without any additional parame- 
ters (dashed line in Fig. ||). However, our data can 
equally well be modelled within SW theory by intro- 
ducing an antiferromagnetic next-nearest-neighbour cou- 
pling of (0.067 ± 0.007) x J, and a priori the experimen- 
tal data cannot be used to discriminate between these 
two possibilities. To help resolve this issue we have un- 
dertaken numerical computations on finite sized systems, 
and do indeed find a ZB dispersion for pure nearest neigh- 
bour coupling |2?j . QMC on lattices up to 32 x 32 spins 
and temperatures down to 0.3 J extrapolate to about 6% 
ZB dispersion p8fl . Exact diagonalisation of an 8 x 8 
system display a ZB dispersion of 4.8%. Both techniques 
yield no ZB dispersion in the 4x4 system, which indicate 
that a significant number of spins are involved in this ef- 
fect. We conclude that the experimentally observed ZB 
dispersion is indeed an intrinsic quantum effect on the 
T = spin dynamics of the pure 2DQHAFSL model with 
nearest- neighbour interactions. Recently a ZB dispersion 
of —13% (i.e. opposite sign) has been found in La2Cu04, 
which is attributed to higher-order spin couplings 29 1. 

Above Tjv, the magnetic excitation spectrum broad- 
ens, making it more difficult to distinguish from the 
phonon scattering. To determine the phonon scatter- 



ing, the SW contribution was simulated by convolving 
the linear SW result with the experimental resolution. 
When subtracted from the 8 K data the remaining inten- 
sity is dominated by phonon scattering. This was scaled 
to the thermal population factor and subtracted from 
the remaining T > T)v data. As this method assumes a 
uniformly renormalised SW dispersion, it obviously does 
not apply to <3||[1,0], and hence only the <3||[1,1] data 
were analysed in this way. Fig. 1(c)- (d) shows examples 
of data sets where the phonon contributions have been 
removed. The background corrected data were analysed 
by fitting cuts along energy to the same parameteriza- 
tion as for the 8 K data, but with the delta function 
replaced by the damped harmonic oscillator (DHO) line 
shape S(u) -uiq) -> ij^r^r^^pr- No clear trends 
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could be found within error for the Q dependence of the 
fitting parameters A, J and T. In the following we there- 
fore choose to discuss the Q averaged values A(T), J(T) 
and f (T). 



T[K] 

Z x 
Z c 
f/J 



8.0 16.2 20.9 25.6 30.4 36.2 45.1 

0.51(4) 0.49(2) 0.51(2) 0.49(2) 0.42(3) 0.58(5) 0.41(8) 

1.18 1.19(1) 1.17(1) 1.17(1) 1.14(2) 1.07(4) 1.03(4) 

0.00(2) 0.08(1) 0.12(1) 0.17(2) 0.24(3) 0.38(5) 0.42(7) 



TABLE I. Temperature dependence of the susceptibility 
and SW velocity renormalization parameters, Z x and Z c , and 
the SW damping F(T)/J. Z c has been normalised to match 
the theoretical value 1.18 at 8 K. The 8 K data were resolution 
limited, and therefore correspond to zero SW damping. 



The amplitude, A(T) = ~(3^z) 2 Z x S, can be used to 
extract the renormalization Z X (T) of the susceptibility. 
(The factor | becomes | for the horizontal banks be- 
low T/v where the moments order almost vertically.) The 
temperature dependence of Z X (T) is listed in Table 1. 
At low T, the experimental value is in perfect agreement 
with the predicted value of 0.51 Q. At higher T, Z x 
appears to decrease, but no clear predictions for the T 
dependence of this quantity have yet been reported. 



The uniform SW renormalization of lvq found along 
[1,1] at 8 K justifies the use of the average J(T), which 
divided by J gives the temperature dependent SW renor- 
malization Z C (T) as shown in Fig. 3. Kaganov and 
Chubukov |^],[30| have calculated the temperature depen- 
dence of the higher-order quantum corrections to SW the- 



ory, obtaining v s (T) — v s (0) (l 



C(3) 
4ir 



where 



£(3) ~ 1.2021. This result is represented by the solid line 
and is in good agreement with our data. The new QMC 
results are indicated by triangles and are also consistent 
with the data. 



3 




0.4 
T/J 

FIG. 3. The SW softening Z c = J/J and damping F/J 
as functions of temperature. The experimental data (cir- 
cles) and new QMC results (triangles) are averages for 
uj > 2 meV. The solid lines are respectively the prediction 
of Kaganov and Chubukov for Z C (T) and the relation 

r(T) = v s (T)/£(T) for the damping. 

Various calculations of the SW damping all predict a 
Q dependence with a minimum at (tt,tt). Around this 
point, the value of T is quite dependent on the defini- 
tion of the line shape. This complication is irrelevant for 
the experimental data, where due to the large incoherent 
background, the fits were restricted to oj > 2 meV ~ J/3. 
The experimentally determined SW damping is shown in 
Fig. ||. QMC results for the damping have been reported 
for 0.35 < T/J < 0.5 We find that these overesti- 
mate the experimental values of T by almost a factor 
of two. To resolve this discrepancy, we have repeated 
the QMC calculations |27j. If the same maximum en- 
tropy (ME) method is used for analytic continuation of 
the imaginary time MC data, then the results of Ref. 
plf are reproduced. We found, however, that it is more 
robust to impose the same DHO line shape as used to 
analyse our experimental data. When this is done the 
QMC values shown in Fig. || are obtained which are in 
good agreement with the experimental data. Scaling ar- 
guments have also been used to calculate the SW damp- 
ing |fl3|] . Since this approach mainly applies to the low- 
energy, long-wavelength behaviour of the system which 
was not probed in our experiments, we believe that it is 
not appropriate to compare it with our data. We do note, 
however, that for T < J/2 the data lie above the predic- 
tion for Q — (7r,7r), and below the prediction extrapo- 
lated to Q = (7r/2,7r/2). Our result for T(T) differ from 
those extracted from NMR in Sr 2 CuCl 2 2 This 
may be becase NMR is a local probe in the low energy 
limit, and the result for T(T) relies on a global asump- 
tion for S(Q,uj). In contrast, our data give a direct and 
asumption free measurement of T(T). 



Spin waves are eigenmodes with respect to antiferro- 
magnetic order and vary on a length scale set by the 
correlation length £(T). The most naive way to estimate 
the lifetime of a SW would be to divide £(T) Q by the 
SW velocity v s (T) fl]|3C]]. In Fig. 4 the solid line is the 
inverse lifetime T(T) = v s (T)/£(T) obtained in this way, 
and is seen to be in surprisingly good agreement with the 
data. 

In summary, we have measured the excitation spec- 
trum of CFTD, which is an excellent physical realization 
of the 2DQHAFSL. Combining our data and numerical 
calculations, the existence of a quantum induced ZB dis- 
persion has been unambiguously established. The finite 
temperature behaviour has been probed up to T > J/2, 
where well-defined SW like excitations persist. The tem- 
perature dependence of the SW softening and damping 
in this strongly quantum mechanical system are remark- 
ably well described, without any adjustable parameters, 
by existing theories based on a quantum renormalization 
of the classical system. 

We gratefully acknowledge Rajiv R. P. Singh for use- 
ful discussions. This work was supported by the Danish 
Research Academy, the UK EPSRC and the EU through 
its TMR and IHP programmes. ORNL is managed for 
the US DOE by UT Battle, LLC, contract DE-AC05- 
00OR22725. 



[1] J. Igarashi, Phys. Rev. B 46, 10763 (1992). 

[2] S. Chakravarty et al, Phys. Rev. B 39, 2344 (1989). 

[3] J.-K. Kim et al, Phys. Rev. Lett. 80, 2705 (1998). 

[4] B. B. Beard et al, Phys. Rev. Lett. 80, 1742 (1998). 

[5] N. Elstner et al, Phys. Rev. Lett. 75, 938 (1995). 

[6] A. Cuccoli et al, Phys. Rev. Lett. 79, 1584 (1997). 

[7] R. J. Birgeneau et al, Phys. Rev. B 59, 13788 (1999). 

[8] M. Greven et al, Z. Phys. B 96, 465 (1995). 

[9] H. M. R0nnow et al, Phys. Rev. Lett. 82, 3152 (1999). 
[10] P. Carretta et al, Phys. Rev. Lett. 84, 366 (2000). 
[11] P. Hasenfratz, Euro. Phys. J. B 13, 11 (2000). 
[12] R. R. P. Singh et al, Phys. Rev. B 52, R15695 (1995). 
[13] S. Tyc et al, Phys. Rev. Lett. 62, 835 (1989). 
[14] A. Auerbach et al, Phys. Rev. Lett. 61, 617 (1988). 
[15] M. Takahashi, Phys. Rev. B 40, 2494 (1989). 
[16] A. Sokol et al, Phys. Rev. Lett. 76, 4416 (1996). 
[17] A. Sherman et al, Phys. Rev. B 60, 10180 (1999). 
[18] S. Winterfeldt et al, Phys. Rev. B 56, 5535 (1997). 
[19] T. Nagao et al, J. Phys. Soc. Jpn. 67, 1029 (1998). 
[20] Y.-J. Wang et al, Phys. Rev. B 56, 10982 (1997). 
[21] M. Makivic et al, Phys. Rev. Lett. 68, 1770 (1992). 
[22] K. R. Thurber et al, Phys. Rev. Lett. 79, 171 (1997). 
[23] S. J. Clarke et al, Solid State Commun. 112, 561 (1999). 
[24] K. Yamagata et al, J. Phys. Soc. Jpn. 50, 421 (1981). 
[25] N. Burger et al, Solid State Commun. 34, 883 (1980). 
[26] Y. J. Kim et al, Phys. Rev. Lett. 83, 852 (1999). 



4 



[27] H. M. R0nnow et ai, To be published. 

[28] O. F. Syljuasen et ai, J. Phys. C 12, LI (2000). 

[29] R. Coldea et ai, cond-mat/006384 (unpublished). 

[30] M. I. Kaganov et al, in Spin waves and magnetic excita- 



tions, ed. A. S. Borovik-Romanov et al. (North-Holland, 
Amsterdam, 1988), Vol. 22.2, p. 1. 
[31] N. Elstner et al, Phys. Rev. B 51, 8984 (1995). 



5 



